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Abstract 

A model for the longitudinal structure function Fl at low x and low Q 2 is presented, 
which includes the kinematical constraint Fl ~ Q 4 as Q 2 — > 0. It is based on the photon- 
gluon fusion mechanism suitably extrapolated to the region of low Q 2 . The contribution of 
quarks having limited transverse momentum is treated phenomenologically assuming that 
it is described by the soft pomeron exchange mechanism. The ratio R = Fl/(F2 — Fl), 
with the F2 appropriately extrapolated to the region of low Q 2 , is also discussed. 



1 On leave from Henryk Niewodniczariski Institute of Nuclear Physics, 31-342 Cracow, Poland. 
2 On leave from the Department of Physics, Jagellonian University, 30-059 Cracow, Poland. 



1 Introduction 



The longitudinal structure function Fl(x, Q 2 ) which corresponds to the interaction of the longi- 
tudinally polarized virtual photon in the one-photon-exchange mechanism of inelastic lepton- 
nucleon scattering is a very interesting dynamical quantity since at least at low x the dominant 
contribution to it comes from gluons. The experimental determination of is in general dif- 
ficult and requires a measurement of the y dependence of the deep inelastic cross section for 
fixed x and Q 2 where, as usual, y = (pip)/(qp),x = Q 2 /(2pq) and Q 2 = —q 2 with pi,p and q 
denoting the four momentum of the incident lepton, the four momentum of the proton and four 
momentum transfer between the leptons respectively, see fig. 1. In the "naive" quark-parton 
model the structure function Fl(x,Q 2 ) vanishes. To be precise, in the limit of large Q 2 it is 
proportional to (< m 2 > + < >)/Q 2 where m is the quark mass and is the transverse 
momentum of the quark which, in the naive parton model, is by definition limited |3[]. 
This remains approximately valid in the leading logarithmic approximation where only collinear 
gluon emissions are considered. The limitation of the quark transverse momentum no longer 
holds in the next-to-leading logarithmic approximation when the point-like QCD interactions 
allow (with probability ~ a s {Q 2 )) < k% > to grow as Q 2 with increasing Q 2 . Thus Fl(x,Q 2 ) 
acquires a leading twist contribution, proportional to a s (Q 2 ) @, |, ||. At small x it is driven by 
the gluons through the g — > qq transition and, in fact, F L can be used as a very useful quantity 
for a direct measurement of the gluon distribution in a nucleon 0. 

In the limit Q 2 — > the structure function Fl has to vanish as Q A (for fixed 2pq). This 
kinematical constraint eliminates potential kinematical singularities at Q 2 = of the hadronic 
tensor W^ v defining the virtual Compton scattering amplitude ||. It also reflects the sim- 
ple physical fact that the total cross section cxl ~ Fl/Q 2 describing the interaction of the 
longitudinally polarized virtual photons has to vanish in the real photoproduction limit. 

The differential cross section for inelastic lepton-nucleon scattering is most frequently ex- 
pressed in terms of the structure function F 2 (x, Q 2 ) and the ratio R(x, Q 2 ) = F L (x, Q 2 )/F T (x, Q 2 ) 
where Ft(x, Q 2 ) = F 2 (x, Q 2 ) — Fl(x, Q 2 ) is the structure function associated with transversely 
polarized virtual photons. Unfortunately measurements of R(x, Q 2 ) are difficult and thus scarce, 
especially at low x and/or low Q 2 . For x < 0.01 and/or Q 2 <0.5 GeV 2 the experimental in- 
formation is particularly poor. The measurements of R are compiled in fig.2. Preliminary 
measurements in the low x region have been announced by the CCFR collaboration |15j and by 



NMC |nj. All the data indicate a small value of R at moderate values of x and an increase of 
R with decreasing x. The data collected in fig.2 come from experiments carried out on different 
targets. This is justified since the measured differences R A — R d do not exhibit any significant 



dependence on x and they are consistent with zero, [p7| . All the data in fig.2 (except those 
of SLAC E140X) have been fitted to a common function (the curve in fig.2), R^ lt , which is 
constructed to extrapolate to theoretically reasonable values outside the kinematical range of 
the data: namely at x — > 0, x — > 1 and Q 2 — > oo. In particular the last constraint means 
that W %t matches R® CD , 0, at high Q 2 . The fit expression should not however be used for 
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Q 2 <0.35GeV 2 which means that no constraints are given on the behaviour of R at low Q 2 . 
This is a serious deficiency in the procedure used to extract polarized and non-polarized struc- 
ture functions from the experimental data, especially in the calculation of radiative corrections. 
Thus the procedure commonly accepted in data analysis is to assume R(x, Q 2 ) = R(x, Q 2 =0.35 
GeV 2 ) with 100% uncertainty for Q 2 <0.35 GeV 2 . The uncertainty of R is then propagated into 
a contribution (usually the dominant one) to the systematic error on the extracted structure 
functions. The need for a more precise determination of R at low Q 2 is thus clearly seen. 

At large Q 2 and large x (x > 0.1) perturbative QCD describes reasonably well the available 
data on the ratio R(x, Q 2 ) JTB|. It has however been noticed that in order to correctly accom- 
modate experimental information from the moderately large Q 2 region, a significant "higher 
twist" contribution, which vanishes as 1/Q 2 at large Q 2 , has to be added as well [ 15] , |H| OT. 



While the longitudinal structure function is (at least theoretically) fairly well understood at 
high Q 2 very little (if anything) is known about its possible extrapolation towards the region 
of low Q 2 . This should be contrasted with the structure function F 2 (x, Q 2 ) for which several 
dynamically motivated extrapolations exist, cf. ||. 

The purpose of this paper is to provide an extrapolation of the structure function Fl(x, Q 2 ) 
towards the region of moderate and low values of Q 2 . This extrapolation will be based on the 
photon-gluon fusion mechanism, essential at low x and suitably extended into the region of low 
Q 2 . A similar description of both the structure functions F2 and Fl has been considered in 
ref. |[21|| . We shall explicitly introduce an additional higher twist term which comes from the 
contribution of the quarks (antiquarks) with limited transverse momentum. We assume that 
this term is described by "soft" pomeron exchange. We treat this contribution phenomenologi- 
cally and determine its normalisation from the (non-perturbative) part of the structure function 
F 2 (x,Q 2 ). The model embodies the kinematical constraint Fl ~ Q A in the limit Q 2 —> (for 
fixed 2pq or v). 

The paper is organised as follows: in the next section we discuss the photon-gluon fusion 
contribution to the structure function Fl and its extension to the low Q 2 region. In section 3 
we present a phenomenological model for the higher twist contribution to Fl at small x which 
is based on the soft pomeron exchange mechanism. Section 4 contains numerical predictions 
for Fl and R, while section 5 contains a brief summary of our results. 



2 The photon - gluon fusion model for the structure 
function Fi 

At small x and at large Q 2 the structure functions can be computed from the kx factorization 
theorem P2L 23, 12JJ] which corresponds to the (virtual) photon - gluon fusion mechanism shown 
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in fig.3. The structure function Fl obtained from the diagram of fig.3 is given by 
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F L (x,Q 2 )=2%Y,e 2 M x >Q 2 ) (!) 



7T 2 ^ q 
q 



where, 

J « (X ' Q2)= J li Jo dP I d2 ^ 2 (l - - ^) 2 /(%, 4) (2) 



and 



D lq = 4 + /3(1 - /3)Q 2 + m 2 , 
D 2q = (k t - fc T ) 2 + 0(1 - /3)Q 2 + m\ . (3) 

a s is the QCD coupling and its argument will be specified later. The transverse momentum 
n' T and the variables y q are defined as follows: 

4 = k t -(1- P)hr (4) 

where Kt is the transverse momentum of the quark(antiquark) in the quark box (see fig.3). 
The variable j3 is the corresponding Sudakov parameter appearing in the quark momentum 
decomposition into the basic lightlike four vectors, p' and q', 

k = x q p' — f3q' + Kt 



X q = X I 1 + 



M 2 x 
q' = q + xp 



P' = P~ -pyTQ (6) 



where M denotes the nucleon mass. The function f(y, k 2 ) is the so - called unintegrated gluon 
distribution. It is related to the (scale dependent) gluon distribution g(y,[i 2 ) by 

V9(V,H*)= r 2 #/(y>4)- (7) 

J f£rp 



The parameters m q are the masses of the quarks. It should be noted that the integrals 
I q defined by (Q) are infrared finite even if we set m q = 0. The non-zero values of the quark 
masses are however necessary if formula (|l|) is extrapolated down to Q 2 = respecting the 
kinematical constraint ~ Q 4 . The non-zero quark masses then play the role of the infrared 
regulator. It should be noticed that the expressions (|],|2]) define Fl as an analytic function 
of Q 2 except for the cut Q 2 < —Am 2 . Equation ([|) can be interpreted as representing the 
rescattering of the qq fluctuation of the virtual photon through the exchange of two and more 
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interacting gluons. Both the diagonal and non-diagonal transitions in this rescattering are 
included. We assume that equations (|I|,0) embody (at least in the average sense and as far as 
the position of the singularities in the Q 2 plane are concerned) also the virtual vector meson 
contributions which couple to virtual photons. This will fix the magnitude of the corresponding 
quark mass parameters to be m 2 ~ m 2 /A where m v denotes the mass of the lightest vector meson 
corresponding to the qq pair. 

The integration limits in (0) are constrained by the condition: 

y q < i (8) 

The condition y q > x is automatically satisfied (see eq.@ ) 

The unintegrated gluon distribution f(x, k 2 ) should in principle be obtained from the solu- 
tion of the Balitzkij-Fadin-Kuraev-Lipatov (BFKL) equation ]26|, |27| . In our calculations we 
shall use the following approximation, 

dg AP (y,Q 2 



dlnQ* 



(9) 



where g AP (y, Q 2 ) satisfies the conventional (LO or NLO) Altarelli-Parisi equations. In this 
approximation one neglects the higher order small x resummation effects in the gluon anomalous 
dimension. The gluon anomalous dimension / -f gg (a s ,u>) is the moment of the splitting function 

j gg (oi s , u) = / dzz"P gg {z, a s ) . (10) 

J 

The BFKL equation effectively resums the leading powers of a s /u i.e. 



)=t+x:A,(t)" (11 > 



w n=A 

(d s = 3a s /ir), where the coefficients D n can be calculated from the kernel of the BFKL equation 
[p8[ . It should be noted that the sum on the right-hand-side of ([TT|) only starts at n — 4. As a 
result it affects the evolution of gluon distribution (and of observable quantities) only at very 
small values of x (x < 10~ 4 ). In our approximation this contribution is neglected. 

Expression ([]]) can be represented in the following "kt factorized" form: 

Fl(x, Q 2 ) =\ X %\ ^fFl(x>, Q 2 , k 2 T )f(^, k 2 T ) (12) 

J X X J Krp X 

where the function F®(x', Q 2 , kj) can be regarded as the longitudinal structure function of the 
off-shell gluon with virtuality k\. It corresponds to the quark box (and crossed-box) diagrams 
in the upper part of the diagram in fig. 3 and is given by the following formula: 

Q 4 V^„2 f 1 jo f i2.J /c ( I f-t , K^+m 2 k?^ 



xa s (3 2 (l-{3) 2 (^--^-) 2 (14) 

where Di q ,D 2q are defined by equation (|3|). The leading twist part of the k t factorization 
formula can be rewritten in a collinear factorized form: 

F L (x, Q 2 ) = / ^C L (x', Q\ a s (Q 2 ))-g(-, Q 2 ) (15) 

Jx X 1 XX 

where the leading powers of a s ln(x'), (ajn(x/x')) have been resummed in Cl(x' ', Q 2 , a s (Q 2 )), 
(x/x'g(x/x' ,Q 2 )). Thus for instance the moment function 

C L {u>, Q 2 , a s {Q 2 )) = f 1 -z"C L (z, Q 2 , a s {Q 2 )) (16) 
Jo z 

is related to the moment function F°(o;,7), by 

C L (u, Q 2 , a s (Q 2 )) = rf gg (a a , u))Ff(u, 7 = l gg (a s , u)) , (17) 

where 



1 /-1/2+ioo (Q 



2\ 7 



n^.Q 2 ) = 7T- , ^(^7) s- as) 

2m Ji/2-ioo V^t/ 



F L V,Q 2 ) = f^F^Q 2 
Jo X 



Our approximation corresponds to setting 7 99 = a s /u. It correctly generates the first three 
corrections to Cl which are presumably most important at moderately small values of x. 

In the "on-shell" approximation one sets k\ = in the argument of F®(x', Q 2 , kj) and 
restricts the integration over k\ to the region k\ <C Kj> ~ Q 2 . After taking into account (|7]) 
the structure function can be expressed in terms of the conventional gluon distribution 

yg(y,Q 2 ), 

I q (x, Q 2 ) = tt £ d(3 J dK ,2 a s (Q 2 )(3 2 (l - (3) 2 ^y q g(y q , Q 2 ) (19) 



where now 



and 



D q = 4 + (3(1 - (3)Q 2 + m 2 (21) 



In order to extrapolate the formula ([!]) to the low Q 2 region (with I q given by (Q) or (|H?D) 
we have to freeze the evolution of g(y, Q 2 ) and freeze the argument of a s (Q 2 ). To this end we 
substitute Q 2 + 4m 2 instead of Q 2 as the argument of a s and of yg(y,Q 2 ). This substitution 
may be justified by analyticity arguments i.e. we want I q to have "threshold" singularities 
for Q 2 < —Am 2 . Moreover, for heavy quarks, it is the heavy mass squared and not Q 2 which 
should control the scale of a s and of the parton distributions for moderately large values of Q 2 . 
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It should be noted that after those modifications the structure function Fl can be continued 
down to Q 2 = 0, respecting the kinematical constraint ( Fl ~ Q 4 in the limit Q 2 — > ). 



It follows from (|T9|) that the structure function Fl in the on-shell approximation can be 
represented by 

mi 



where 



0:.s 



dy ( x^ 



x 



3// \ Q 2 (?/-a;) 



4m ^ y*?(y> Q 2 ) 



O'.s 



dy ( x" 



( 



In 



and 



7T ^, y vy, 



+ v /l- 


Arri^x \ 

Q 2 (y-z) 








Am 2 








Q 2) ' 





yg(y,Q' 



(22) 



(23) 



(24) 



(25) 



It should be noted that in the large Q 2 limit, the first term becomes equal to the standard 
QCD formula describing the leading order contribution to Fl originating from the photon-gluon 
fusion (§§, 

F L (x, Q 2 ) = 2^ £ e 2 / 1 (-) 2 ( 1 - -)V9(V, Q 2 ) ■ (26) 
Y q Jx y y 

In the large Q 2 region formula (|22|) with non-zero quark masses generates also power corrections 
enhanced by logarithmic factor present in the term jj- 2 \ 

When analysing the kx factorization formula it is convenient to use the "on-shell" approx- 
imation in the region k 2 < &q, where k^ is a parameter of the order 1 GeV 2 and to leave the 
remaining part (k 2 > fcg) unchanged. The first term is then given by (0), (f23"l), (^4j) where 
gluon distribution is evaluated at the scale k^ instead of Q 2 p4]| . This gives the following 
representation for the structure function Fl'- 



F L (x,Q 2 ) = F L0 (x,Q 2 ;k 2 T0 ) + 



2|fE^/ ffe(4-4o) [d? J d 2 K' T a s (3 2 (l-P) 2l -( 



I, 1 



D 



D 



2q 



2 f(y q ,k. 



T) 



(27) 



where Flo(x, Q 2 ; k T0 ) is given by equations (|22"D - (p4]) with 7/g(|/, /cqt) instead of yg(y, Q 2 ) in 
the integrals ( p3|) and fl2"4"|). The argument of a s in the formulae (0), ( p4|) defining Ft, (a?, Q 2 ) = 
Flo{x, Q 2 ', k TQ ) and in the second term of Q27D will be set equal to Q 2 + Am 2 . 
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3 A model for the higher twist contribution to Fi 



The phenomenological description of the experimental data for the ratio R(x, Q 2 ) in the region 
of moderately large values of Q 2 and large values of x (x > 0.1) requires a significant higher 
twist contribution, i.e. the term which vanishes as 1/Q 2 for Q 2 



oo 



In this section we present the phenomenological model for Ft which takes into account the 
higher twist effect. The idea is to treat the contributions from the regions of low and high 
values of quark transverse momenta in different ways. That is we divide the integration over 
Krp into two parts: the region < < k' 2 t , and k'^ > k'q T , where K® T is an arbitrary cut-off, 
chosen to be of the order of 1 GeV 2 . 



We leave unchanged the contribution coming from the high k'^ region whereas in the low n'f 1 
region, which is presumably dominated by the soft physics, we use the "on-shell" approximation 
and make the substitution: 

a s zg(z,Q 2 ) -> A . (28) 
This gives the following representation of the higher twist contribution to Ft: 



2^e^ f'd^il-P) 2 f" 

„ 7T JO JO 



,12 

"^1 



(29) 



where D q has been defined in ([H]). The constant A is not a free parameter. We estimate it 
from F 2 assuming that the non-perturbative contribution to F 2 in the scaling region also comes 
from the region of low values of k,'^ and is controlled by the same parameter. The term F^ T 
can be interpreted as representing the contribution of soft pomeron exchange with intercept 
equal to 1. 



It should be noted that F^ T given by equation (|29|) will vanish as 1/Q 2 in the high Q 2 limit 
(modulo logarithmically varying factors). We call it therefore a "higher twist contribution". It 
should be noticed that this term will also respect the kinematical constraint Ft ~ Q 4 in the 
limit Q 2 -> 0. 



The corresponding formula describing the photon 
function F T is 



gluon contribution to the structure 



F T (x, Q 2 ) 



r, \ "* 2 Q 2 

q 



/ dp dn' 2 - f g(- f ,Q 
Jo J x' x' 

p 2 + (1 - P) 2 



X 



2 K>rr 



2 thrpfcrp 



+ 



(30) 



We integrate this expression over the low region (n® < k'q T ) after substituting a s zg — > A and 
retain its limiting value as Q 2 — > oo. We identify this contribution to Ft with the "background" 



term F 2 Bg of ref.Eg 



i.e.: 



F? 9 



T e 2 

7T 



dt 



.12 
"OT 




m 2 K'^ 
n 



(31) 



where now 



m 2 + n' 2 + t . 



(32) 



A possible weak x dependence of F^ g is neglected and we set F 2 B9 = 0.4. 



The complete structure function Fl is represented as Fl = F^ T + F[ T where F^ T is 
calculated from equations (|22|) - (|24]) , or from equation (|27|) with the additional constraint 
Kj- > k'q T in order to avoid double counting. 



4 Numerical results 

In this section we present the numerical analysis of the structure function Fl- We shall also 
perform the analysis of the ratio R(x, Q 2 ) and compare it with the available data. 

In our calculations we have included the contributions of the u, d, s, c quarks with masses 
0.35, 0.35, 0.5, 1.5Ge^ respectively. The parameter k' 2 t was set to be equal 0.8 GeV 2 . We have 
varied this parameter in the interval 0.8-1.5 GeV 2 and found that the results for R are not very 
sensitive to these changes, i.e. the ratio R does not change by more than 10%. 

We have used the GRV [0 and MRS (A) @ parton distributions both in the "on-shell" 
and "off-shell" approximation. In the case of the GRV parametrisation LO gluons and QCD 
coupling constant were used. In our calculations we froze the evolution in Q 2 in the argument 
of a s (Q 2 ) and also in the function yg(y,Q 2 ) as explained before. Besides the dominant (at 
small x) photon-gluon fusion — > qq contribution, we have also included in Fl a contribution 
from quarks taking into account threshold effects in the corresponding convolution integral i.e. 

AF L (x, Q 2 ) = £ 2 <*.(Q 2 + 4™;) 4 f 1 *L(Z.)* [<h (y, Q 2 + 4m 2 ) + Q 2 + 4m 2 )] . (33) 

where x q is given by (p5|). 

In figures 4 and 5 we show our results for the longitudinal structure function calculated 
without an additional higher twist contribution. Fl is plotted as the function of Q 2 for different 
values of x, and as a function of x for different values of Q 2 . The theoretical curves reflect the 
different behaviour of the MRS (A) and GRV gluon distributions. The difference between the 
on-shell and off-shell approaches is in general small, particularly in the small Q 2 region. It can 
be seen from fig. 5, that the magnitude of Fl strongly decreases as Q 2 — > 0. 

In figures 6 and 7 we show results which include the additional higher twist contribution. It 
can be seen from these plots that the resulting magnitude of the longitudinal structure function 
is now bigger in the low Q 2 region than in the previous case. It should also be noted that the 
structure function Fl flattens as a function of x in the low x and low Q 2 region. This is a direct 
consequence of the "soft pomeron" parametrization of the higher twist term. 
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We have also calculated the ratio R(x, Q 2 ) = F L (x, Q 2 ) / (F 2 (x, Q 2 ) — F L (x, Q 2 ))- The struc- 
ture function F 2 was parametrized following the model of ref . . In this model the structure 
function F 2 is represented by 

+ + (34) 

The first term describes the VMD contribution and the sum extends over the low mass vector 
mesons. The parameter 7„ can be estimated from the leptonic widths of the vector mesons, 
a v (s) is the total vector meson - proton cross section, and s = 2pq — Q 2 + M 2 where M is the 
nucleon mass. The structure function F 2 (x,Q 2 ) is obtained from the QCD improved parton 
model. The variable x is defined as 



x = x(l + |j) . (35) 



In the case of the GRV parametrization only the light quarks contribution was included in 
F 2 S in the formula (J33J). The charm quark contribution to F 2 was calculated from the on-shell 
photon-gluon fusion. 

In figures 8 and 9 we show our results for the ratio R(x,Q 2 ). We see that this quantity is 
sensitive to the different parton distributions which were used in calculating both Fl and F 2 . In 
particular the structure function F 2 obtained from the formula (|3~4]) with F 2 AS calculated from 
the MRS (A) parametrization turns out to be steeper function of x in the low Q 2 region than 
that when F 2 S is calculated from the GRV partons. As a result the ratio R decreases with x as 
x — > 0. This effect is weaker if Fl does not contain the higher twist modification. The spread 
of the theoretical predictions can be used as an estimate of the corresponding errors caused by 
the uncertainty of R in the low Q 2 , low x region. For the intermediate values of Q 2 this model 
coincides with SLAC parametrization, whereas for lower values of Q 2 , outside the range of the 
validity of this parametrization, it quickly vanishes as Q 2 — > 0. 

In fig. 10 we show extrapolation of our model calculations for R up to x = 0.1 and confront 
this extrapolation with the available data. Although this value of x is already too large for 
our model to be applicable one can see that it predicts reasonably well the magnitude and the 
Q 2 dependence of R. The agreement with experimental results improves if we renormalise the 
phenomenological parameter A by a factor of 1.2 - 1.5 . 



5 Summary and conclusions 

In this paper we have presented possible parametrization of the longitudinal structure function 
Fl and of the ratio R(x, Q 2 ) in the region of low Q 2 and low x. We based this parametrization 
on the photon-gluon fusion mechanism suitably extrapolated into the region of low Q 2 . The 
extrapolation respected the kinematical constraint Fl ~ Q 4 in the limit Q 2 — > 0. We have also 
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included a separate contribution from the higher twist term which corresponds to soft pomeron 
exchange. Its coupling to external virtual photons was assumed to be proportional to that of the 
(on-shell) gluons but with limited transverse momentum of the quarks (antiquarks) within the 
quark box (and crossed-box) diagram, eq. (EH). We have also calculated the ratio R(x,Q 2 ) = 
Fl(x, Q 2 )/Ft(x, Q 2 ) utilising the parametrization of the structure function F 2 (x, Q 2 ) from ref. 
p2[ . Our results show that the structure function Fl and the ratio R become negligibly small 
in the region of low Q 2 (Q 2 <0.1 GeV 2 or so). The ratio R(x,Q 2 ) turns out to be essentially 
independent of x in the low x, low Q 2 region. We have also noticed that in this region the 
magnitude of both Fl and of the ratio R is relatively insensitive to the variations of the input 
parton distributions. We believe that our results will also be useful in the radiative corrections 
procedure as well as in extracting the magnitude of the spin dependent structure functions 
from the experimental data. The codes for calculating R(x, Q 2 ) are available upon request 
from a.m.stasto@durham.ac.uk. 
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Figure captions 



1. Kinematics of inelastic charged lepton-proton scattering in the one-photon-exchange 
approximation and its relation through the optical theorem to Compton scattering for 
the virtual photon; pi, p and q denote the four momenta of the incident lepton, proton 
and virtual photon respectively. 

2. Compilation of measurements of R; data come from SLAC E140X ||, SLAC E140 and 
the global analysis of earlier SLAC experiments P|, CDHSW pf, BCDMS ITJ and 



EMC |I3| experiments. Errors show the statistical and systematic uncertainties added in 
quadrature but in case of SLAC E140X and the global analysis of the SLAC data they 
do not include the systematic error due to radiative corrections (equal to approximately 
0.012 and 0.025 respectively). The curve is a fit to all of the data (excluding E140X), 



[ 1~0| , jTJJ. The figure is taken from [ |15|j . 



3. The diagramatic representation of the photon-gluon fusion mechanism and of the 
factorization formula (IT! 



4. Longitudinal structure function Fi calculated as a function of x for different values of 
Q 2 assuming the MRS (A) and GRV gluon distributions for both the on- and off-shell 
approaches. 

5. As fig.4 but F L is calculated as a function of Q 2 for different values of x. 

6. As fig.4 but results include a contribution from higher twist. 

7. As fig. 5 but results include a contribution from higher twist. 
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8. R(x, Q 2 ) plotted as a function of x for different values of Q 2 and containing a contribution 
from higher twist. 

9. As fig. 8 but R is plotted as a function of Q 2 for different values of x. 

10. Comparison of the model calculation for R, extrapolated to x—0.1, with the data. The 
higher twist is taken into account in the calculations. See fig. 2 for a description of the 
data. 
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